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Microscopic mechanisms of the puzzling insulating ferromagnetism of half- filled La4Ba2Cu20io are 
elucidated with energy-resolved Wannier states. The dominant magnetic coupling, revealed through 
evaluated parameters (£, U, and J), turns out to be the intersite direct exchange, a currently ignored 
mechanism that overwhelms the antiferromagnetic superexchange. By contrast, the isostructural 
Nd4Ba2Cu2 0io develops the observed antiferromagnetic order via its characteristics of a ID chain. 
Surprisingly, the in-plane order of both cases is not controlled by coupling between nearest neighbors. 
An intriguing pressure- induced ferromagnetic to antiferromagnetic transition is predicted. 



Unlike most half-filled cuprates that feature the anti- 
ferromagnetic (AF) ground state as expected from the 
Hubbard model, few of the cuprates possess exceptional 
ferromagnetic (FM) order instead ||, [|]. In par- 
ticular, half-filled Cu-d orbitals in the "brown phase" 
La4Ba2Cu2 0io (La422) is found to align ferromagneti- 
cally below 5K ^| , in great contrast to the isostructural 
(Fig. |l|) Nd 4 Ba 2 Cu20io (Nd422) that retains antiferro- 
magnetism below 7.5K ||, ||, ||. Such an unusual insulat- 
ing FM phase poses a great challenge to our quantitative 
theoretical understanding of microscopic mechanisms in- 
volved in real materials. 

Up to now, there are only limited pioneering attempts 
at identifying the quantum processes responsible for this 
intriguing behavior of La422. With elaborate perturba- 
tion on a multiband Hubbard-like model 0, it is sug- 
gested that destructive interference of hopping paths (be- 
tween 6 th order and 8 th order terms) may suppress AF 
coupling and give a small FM coupling along the z-axis in 
a very narrow parameter range, while the Goodenough 
process [|| between nearest neighbors produces the in- 
plane FM order. (More recently, the possible "crude 
link" to the "flat-band ferromagnetism" from the Hub- 
bard model was also pointed out ||.) However, this de- 
scription is not very satisfactory; besides the problem of 
the narrow parameter range, the resulting strength of the 
FM coupling is too small to account for the experiments 
[Q. Furthermore, the assumption of FM coupling be- 
tween the nearest neighbors in the plane is apparently in 
contradiction to the AF order in Nd422 ||. 

From the theoretical point of view, descriptions of 
quantum magnetism have been given mainly via phe- 
nomenological models like the Stoner, Heisenberg, Hub- 
bard, or other tight-binding models (8[ [llj, with ad- 
justable parameters to fit the experiments. These ap- 
proaches, though intuitive and computationally manage- 
able, lack the quantitative detail of the complex interac- 
tions that occur in real materials. Sometimes even the 
sign of spin-spin coupling is simply chosen to fit the ex- 
perimental magnetic order. On the other hand, a normal 
practice of more quantitative ah initio approaches [ p"T| , 
involving comparing total energy corresponding to differ- 



ent magnetic orders within the density functional theory 
(DFT) |l2j , not only suffers from the uncertain quality of 
existing approximate energy functionals jl2| for the very 
sensitive magnetic systems, but also hides all the micro- 
scopic processes in a black box. The development of a 
unified scheme that simultaneously gives intuitive micro- 
scopic insight and quantitative realism of ah initio level 
is thus highly desired. 

In this letter, aiming to resolve specifically the micro- 
scopic mechanisms of the unusual insulating ferromag- 
netism in La422, we attempt the first step of one such 
ambitious scheme by evaluating parameters (t, J7, and J) 
of (reformulated) ah initio second quantized Hamiltonian 
in the basis of localized energy-resolved Wannier States 
(WSs) |l3|, |l|, JTl [L6|, [L7|, constructed from all-electron 
DFT orbitals [O. A clear picture of microscopic pro- 
cesses involved naturally emerges from these parameters 
and the spatial distribution of the WSs. The dominant 
mechanism turns out to be the intersite FM "direct ex- 
change" jlO| that is currently ignored in the microscopic 
studies of this system. This process, occurring mainly at 
the La and O sites, overwhelms the weak tendency to- 
ward AF order via superexchange of Hubbard type . 
The same analysis is then applied to Nd422 and used to 
demonstrate its characteristics of ID AF chain, as ex- 
perimentally observed ||, ||. Surprisingly, in both com- 
pounds, the dominant (FM) exchange is found not with 
the nearest neighbors (currently assumed), but with sites 
above/below them, which simultaneously generates the 
proper in-plane order. Finally, the crucial role of "chem- 
ical" effect is illustrated with a numerical simulation of 
La422 "under pressure", which suggests an intriguing 
pressure-induced FM to AF transition. Predictions are 
made for further experimental investigation. 

In order to uncover the underlying many-body interac- 
tions responsible for the observed magnetic order, map- 
ping the system onto a simplified effective model Hamil- 
tonian (Tt], [l^] is avoided. Instead, attention is paid to 
the full ah initio Hamiltonian, rigorously reformulated 
into a "fluctuation" form that explicitly utilizes DFT so- 
lutions p9 (without "double-counting" many-body ef- 
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fects): 

-f other (c^Cjj'a) {c&p c D'j^j terms + constant terms, (1) 

where summation over barred variables is understood and 
= 0— {o) DFT denotes "fluctuation" from expectation 
value (o) DFT of operator o within DFT. Notice that the 
direct exchange (the 5 th term) that directly couples spins 
is omitted in all the existing Hubbard-model based anal- 
ysis of La422, and turns out to play a crucial role. 

Extraction of parameters of this full Hamiltonian via 
traditional constraint DFT calculations [[l7[ |l9| is ex- 
tremely complicated (if even possible.) Instead, "bare" 
parameters t^ v = — (/i \h DFT — v xc \ v) (1 — 8^) + 

J,Act^) DFT ~ -( M |^ FT |^)(i-V). V = 

(fiu \v\ v/j) (1 — and = (fifi \v \ fifi) are directly 
evaluated on the basis of half-filled, low-energy, local- 
ized WSs that possess the spin moment. (h DFT ', v xc , 
and v are the DFT Hamiltonian, exchange-correlation 
potential, and bare Coulomb interaction.) These param- 
eters prove to be insightful to illustrate the insulating 
ferromagnetism of interest, as they control how electrons 
(and spins) interact with each other (in a many-body 
environment) in real materials. Now, since solving the 
multiband Eq. (Q) currently remains an important and 
yet challenging task, we proceed by identifying two pri- 
mary microscopic mechanisms that appear naturally in 
Eq. ([j]): Hubbard- type AF superexchange Q and FM 
direct exchange [|l(| (between the WSs). Their strength 
can be estimated by summing phenomenological formu- 
lae: J sx r 4t?./Wj and J DX ~ J^/2, respectively, 

over leading exchange paths, where i and j are site index 
and Wj = Uj/2 approximately accounts for effects of on- 
site repulsion, competing off-site Coulomb repulsion, and 
other possible virtual processes involving states of higher 
energy (screening.) (The factor in Wj is chosen to give 
roughly consistent size of the expected Mott-Hubbard 
gap in cuprates. Considering the intersite nature of the 
direct exchange, the factor in J DX is reasonable [ |30| . 
More accurate estimation of these factors will not qual- 
itatively change our conclusion on the insulating ferro- 
magnetism of La422.) 

The employment of the WSs (or "molecular orbitals 
in crystal") [lj, ||, as basis, instead of other 
orthonormal complete sets, is strongly physically moti- 
vated. In addition to full awareness of crystal and local 
symmetries j JT5| , the constructed WSs are the most local- 
ized within the sub space of low- energy excitation, which 
facilitate direct physical interpretation consistent with 
localized spins interacting with each other. This "en- 
ergy resolution" naturally follows the identity of the sub- 
space spanned by our WSs and that by the chosen eigen- 
states of h DFT . In addition, single-particle properties of 



these WSs — ( "self-interaction" |2l| free) orbital energy 
£fia = (ft \h DFT - v xc \ /i) - Ufj, (n m ) DFT , occupation 
number, and spin moment — are directly translated with 
information from the DFT calculation, and are known 
before solving the many-body Hamiltonian. Specifically 
for La422, we are thus able to focus only on the half- filled 
low-energy WSs corresponding to the two partially occu- 
pied bands across the Fermi energy, as other particle/hole 
states are of much higher energy and with zero spin mo- 
ment. Furthermore, superexchange involving other states 
diminishes due to orthogonality between subspaces. 

In this work, h DFT is chosen via LDA [|| for La422 
and via LSDA+U || for Nd422 (to prevent Nd /-states 
from incorrectly falling on the Fermi energy.) Values of 
U and J (chosen to be 6.7 and 0.7 eV) in LSDA+U, 
employed to split the /-states, are not crucial as long as 
their hybridization with the WSs is eliminated. 

The WSs in unit cell R and of orbital m are con- 
structed by \Rm) = \km) e~ lk ' R / \/N , where N is the 
number of discrete /c-points in the reciprocal space, and 
| km) = \ipkm') (ipkfh' I fern) (band- mixed Bloch state of 
crystal momentum k) is transformed from the all-electron 
[ p4| DFT eigenstate \^km') [^5j of band index m' , to 
facilitate maximum localization of the WSs. We fol- 
low the first step in Ref. |l4| to calculate (tpkm'lkm) by 
projecting onto the chosen subspace of \ijjkm') (belong- 
ing to the two partially occupied bands) narrow Gaus- 
sian states, \g m ), each centered at one of the main lobes 
of Cu-d contribution, followed by symmetric orthogo- 
nalizaton: {^km^km) = (i/>km'\9m") Mfh» m ] M m"m 
(gm"\ipkfh') (i/Jkm'\9m) • Due to the concentrated con- 
tribution from Cu-d orbitals and the fact that there is 
only one low-energy WS per Cu site, the resulting WSs 
are practically the same as the maximally localized WSs 
constructed with the complete iteration procedure 
While denser /c-sampling can slightly refine the shape of 
the resulting WSs, a 7 x 7 x 8 lattice is good enough 
to give reliable description, as the "aliasing" caused by 
overlap of WS with itself due to periodicity is eliminated. 

The resulting low-energy WSs of La422, each centered 
at one Cu site, are shown in Fig. |l[ in which isosurfaces 
of \(x\Rm)\ 2 in real space x are plotted, along with the 
crystal structure. Notice that the "hybridization" be- 
tween Cu-d, O-p and La-d orbitals is naturally built into 
each WS; in the traditional atomic-orbital based anal- 
ysis |l7| , such hybridization would require inclusion of 
complicated hopping and interaction between several or- 
bitals from Cu, O, and La covering larger energy, which in 
turn obscures straightforward visualization of key physi- 
cal processes. The anti-bonding nature of the Cu-d and 
O-p orbitals is also easily observed, as well as the fasci- 
nating "ring" structure at the La sites, which turns out 
to play an important role. 

The spatial distribution of the WSs also reflects that 
of the spin density: (x\Rm) (4^%^/ a ) {R'm'\x) - 
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FIG. 1: Crystal structure of La422 (a) and its top view (d) 
shown within 2x2x2 unit cells. Most atoms outside the 
1 st unit cell are removed from panel a) for clarity, b) and c): 
illustration of low-energy Wannier states with R = (0, 0, 0). 

| (x\Rfh) | 2 (n^rna) for localized WSs. We found that 
about 50% of the spin-moment resides at the Cu site, 
while 10% at each of the four O near Cu and 2% (x2 
from two WSs) at each La; the last naturally accounts 
for the experimentally observed hyperfine field for La [|| . 
Our WSs also illustrate, for the first time, the shape of 
the spin density at the La sites (as superposition of d 
orbitals of various symmetry), which has not been com- 
pletely resolved by NMR measurement 0. 

Figure || lists our calculated leading parameters. Note 
that tR m ,R'm' ~ (Rm\km ff ) (kfh ff \R f m f ) is easily 
evaluated with DFT eigenenergy Skm" , in contrast to the 
conventional procedure of tight-binding fit |26|. (Since 
the resulting Uj corresponds exactly to the original DFT 
eigenenergy, this can be viewed as perfect " downfolding" 
p7|.) Ui, and Jij are obtained via careful numerical inte- 



gration [|20|, |28[] , instead of parameterized Slater integrals 
p9[ that assume atomic orbitals with well-defined an- 
gular momentum. The resulting Ui ~ 7.5 eV is much 
smaller than that obtained from atomic Cu-d orbitals, 
reflecting non-negligible hybridization. 

The microscopic mechanisms for the unusual FM or- 
der are now easily visualized. Consistent with the ex- 
tremely narrow bandwidth of the low energy bands, hop- 
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FIG. 2: Calculated parameters (in meV) for La422 and 
schematics of leading hopping (blue single arrows) and direct 
exchange (red double arrows) channels. Coordinate system 
for site indices i and j are defined for convenience. 



ping (tij) between the WSs is very small (Fig. ||), of 
a few meV, which brings about a weak AF tendency 
via superexchange ]l8p . Along z-axis (1 st table), the 
FM direct exchange, commonly discarded in Hubbard- 
like models, clearly overwhelms the AF superexchange 
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This process occurs mainly at the 



"ring" structure of WSs at the La sites (Fig. Q), where 
large hyperfine field occurs Q . In the language of atomic 
physics, this direct exchange is similar to what gives the 
first Hund's rule at La sites. 

Interestingly, instead of the nearest neighbor interac- 
tion currently assumed by existing studies, the in-plane 
FM coupling (2 nd table) turns out to be dominated by 
the (101) direct exchange {above / below the nearest neigh- 
bors), occurring mainly at O sites near Cu (Fig. [l[) Hav- 
ing the interplane FM alignment, a strong in-plane FM 
coupling is thus established, unchallenged by any AF cou- 
pling, due to zero hopping dictated by the local symmetry 
[ |l5[ . The resulting coupling constants are of the same or- 
der as the observed transition temperature (~5K), which 
further rationalizes our phenomenological estimation. 

For comparison, the AF order in Nd422 is also analyzed 
with the same technique. (According to the analysis of 
neutron diffraction measurement |^|, the Nd moment is 
perpendicular to and "disconnected" from Cu moment, 
and thus is excluded in our analysis.) One immediately 
notices in Fig. || that the (001) hopping is considerably 
(~10 times) larger, which in turn greatly enhances the 
AF superexchange and thus allows it to overcome the 
(slightly enhanced) FM direct coupling along the z-axis. 
Such strong ID characteristics is in excellent agreement 
with the optical and neutron measurement [|[ ||. In the 
xy-plane, however, there exits no AF coupling between 
nearest neighbors but only FM exchange, primarily along 
(101), as in the case of La422. (The competing (100) FM 
coupling is weaker, and connects to fewer neighbors.) We 
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FIG. 3: Upper panel: calculated parameters (in meV) with 
low-energy WSs of Nd422 (c.f.: Fig. ^|). Lower panel: com- 
parison of WSs of Nd422, La422, and La422 under pressure. 
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thus reach a surprising conclusion that, having the AF 
order between planes, the observed in-plane AF order is 
actually caused by the (101) FM coupling. 

Insight into the dramatic difference in the z-axis hop- 
ping in these two compounds can be obtained by compar- 
ing the corresponding WSs, shown in Fig. [|. In Nd422, 
the development of extra p-contribution (circled) at O 
sites that belong to the WS one layer below/above, pro- 
vides additional efficient hopping channels along the z- 
axis. By contrast, the weak z-axis hopping in La422 is 
mainly through the (inefficient) "ring" structure at the 
La site, and is further weakened by the residual hop- 
ping via O sites (not shown in Fig. ||), due to the phase 
difference between these two different routes. (This dif- 
ference in hopping channels is clearly reflected in the sign 
of (001) hopping, since Uj, unlike J^, is sensitive to the 
phase shift.) This weakening effect resembles the "inter- 
ference" discussed in Ref. [Q , except that in our analysis 
it is not responsible for the observed FM order. 

The driving mechanism of this spatial redistribution 
of charge/spin is numerically examined by another set of 
calculations for La422 with lattice constants and atomic 
positions identical to Nd422 (La422 "under pressure" 
with ~ 2% reduction of the lattice constants.) The re- 
sulting WSs, shown in Fig. ||, do not possess the ex- 
tra O-p structure. Consistently, the key z-axis hopping 
only reaches 30% of that in Nd422 and the sign of (001) 
hopping remains positive. Evidently, the charge redis- 
tribution is mainly driven by the change of chemical en- 
vironment surrounding Nd, which is smaller in size and 
consists of 3 /-electrons and 3 more protons. 

Intriguingly, while the resulting magnetic order re- 
mains FM, in good agreement with experiment on La422 
under pressure || , moderate enhancement of Uj and Jij 
is observed due to slightly larger covalency, i.e. overlap 
between WSs. As a result, AF superexchange starts to 
approach the strength of the FM direct exchange (since 
t is squared) indicating a FM to AF transition upon 
higher pressure. We propose further experiment with 
pressure slightly higher than 9GPa to observe this fas- 
cinating FM to AF transition. Alternatively, this tran- 
sition should be more easily achievable with moderate 
pressure on (Lao.sNdo.2)4Ba2Cu20io. 

In summary, microscopic mechanisms of the puzzling 
insulating FM order observed in La422 are identified 
based on calculated parameters of interactions between 
half-filled localized WSs. The long-omitted FM direct ex- 
change is shown to overwhelm AF superexchange. The 
spatial distribution of the constructed low-energy all- 
electron WSs provides detailed insight into the exchange 
processes and the partial contribution of spin- moment, 
in agreement with the NMR measurement. The isostruc- 
tural Nd422 is shown to possess characteristics of ID 
AF chain along z-axis, as experimentally observed, via 
extra efficient hopping through O sites. The in-plane 
magnetic order of both compounds turns out to be intro- 



duced via dominant FM coupling with (101) neighbors, 
not the nearest ones currently assumed. The crucial role 
of chemical replacement is illustrated via numerical test 
on La422 "under pressure" , which further predicts an in- 
triguing pressure-induced FM to AF transition. While 
future development for the many-body solution is de- 
sired, our scheme already allows detailed, quantitative, 
and intuitive description of quantum magnetism that 
complements the common model approaches. 
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